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Abstract: A general one-loop scattering amplitude may be expanded in terms of master integrals. The 
t — coefficients of the master integrals can be obtained from tree-level input in a two-step process. First, use 



known formulas to write the coefficients of (4 — 2e)-dimensional master integrals; these formulas depend on 



an additional variable, u, which encodes the dimensional shift. Second, convert the u-dependent coefficients 
of (4 — 2e)-dimensional master integrals to explicit coefficients of dimensionally shifted master integrals. 
This procedure requires the initial formulas for coefficients to have polynomial dependence on u. Here, 
we give a proof of this property in the case of massless propagators. The proof is constructive. Thus, as 
a byproduct, we produce different algebraic expressions for the scalar integral coefficients, in which the 
polynomial property is apparent. In these formulas, the box and pentagon contributions are separated 
. , explicitly. 
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1. Introduction 



Detailed calculations of multi-particle scattering events are needed in order to analyze new physics 
at the experiments of the Large Hadron Collider. Computational complexity increases rapidly with the 
number of legs, even at the amplitude level. New and improved algorithms are being developed to meet 
these needs. Recent progress at next-to- leading order has been reviewed in 

Scattering amplitudes at one-loop level can be understood in terms of an expansion in master integrals 
[§, HJ. The coefficients of the master integrals may be obtained by direct reduction, or alternatively by 
solving constraint equations derived from singular structures, most notably unitarity cuts |4], ||, ||, [5], ||, ||, [K| 
11|, 0, g]| 0, H, 0, 0, Hl> Hl> S H [2|, |§ g|, |§ |27|, |28], ||, |(| . In order to obtain complete physical 
amplitudes from unitarity cuts, we can work in dimensional regularization, where D = 4— 2e [31, [32], [3^, 34]. 



By now, explicit analytic formulas for these coefficients are available [20, 23, 27, 28]. The input quantities 
are taken from the complete tree- level amplitudes involved in unitarity cuts. There are other promising 
algorithms for finding the coefficients in 4 or D dimensions |3^, 17, 24, EjJ, or specifically the additional 



"rational" parts supplementing a pure 4-dimensional expansion [p6|, 37 , 38 



41, 42, 43, 44, 30 



The formulas of [27|, developed in the context of the .D-dimensional unitarity analysis of [18, 21 1, are 
coefficients of (4 — 2e)-dimensional master integrals; these formulas depend on an additional variable, u, 
which encodes the dimensional shift. To finish the calculation, we convert the u-dependent coefficients of 
(4 — 2e)-dimensional master integrals to explicit coefficients of dimensionally shifted master integrals. 1 

We are presently concerned with the adaptation of the formulas of |27j to an efficient numerical 
algorithm. Two particular issues are addressed in this paper: 

• Because the coefficients of the (4 — 2e)-dimensional integrals are polynomials in the variable u, a 
direct numerical implementation is not obvious. 

• The algebraic expression of boxes includes both box and pentagon contributions. The pentagon 
contribution is signaled by the {au + b) factor in the denominator. 

Our aim is to solve these two problems. More concretely, in this paper we accomplish the following: 



x As an alternative to this last step, complete coefficients of (4 — 2e)-dimensional master integrals could be obtained with 
the recursion and reduction formulas of fli] . BlJ . 
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The proof of the polynomial property of u: In previous work, some evidence for this assumption was 
provided. Now, we give a complete proof. 



• Simplifying our previous expressions: The algebraic expressions for coefficients given in [£7| were the 
full polynomials in u, i.e. a sum of terms of the form c n u n . Here, we give expressions for evaluating 
c n directly from input quantities. 

• Separating the coefficients of boxes and pentagons: We give explicit, separate expressions for coeffi- 
cients of boxes and pentagons. 

For simplicity, the results here are specific to amplitudes with massless propagators. Generalization 
to the massive case is straightforward for the coefficients of master integrals that have nonvanishing cuts. 



Based on the present paper, the generalization to the massive case has been presented in [46]. We work 



within the spinor formalism [47], reviewed in |4J 

The paper is organized as follows. In Section 2, we organize our input quantities from tree amplitudes, 
define some key vectors and spinors from the input, and briefly discuss the dimensional shift. Then we 
proceed to the simplifications of the formulas for coefficients, and the proofs that they are polynomials in 
u. These are given in Sections 3 and 4 for triangles and bubbles, respectively. In Section 5, we address box 
coefficients, and for the first time we present separate formulas for box and pentagon coefficients. Section 
6 contains an application of these formulas, within the example of the 5-gluon amplitude. In Section 7, 
we close with a discussion and comparison to a couple of other recent approaches to the problem of one- 
loop amplitudes. Appendix A contains our definitions of master integrals and dimensional shift identities. 
Appendix B contains alternate, more explicit expressions for the triangle coefficients which may be better 
suited for numerical evaluation, since the derivatives have been taken analytically in every case that arises 
in a renormalizable theory. Appendix C contains many of the details of the polynomial proof for bubble 
coefficients. Appendix D contains analytic expressions used in cuts of pentagons. 



2. Setup and Definitions 

In this section, we set up some key conventions and definitions used in expressing the coefficients of master 
integrals, and in our proofs of polynomial dependence. 

2.1 Unitarity method 

The unitarity cut of a one-loop amplitude is its discontinuity across a branch cut in a kinematic region 
selecting a particular momentum channel. Specifically, we denote the momentum vector by K. Then, K 2 
should be positive, and all other momentum invariants should be negative. The vector K will be a sum of 
momenta of some of the external legs. The discontinuity is given by 
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where the Lorentz-invariant phase-space (LIPS) of a double cut is defined by inserting two <5-functions 
representing the cut conditions: 2 



d D <5> = d D p 5(p 2 )5({p - Kf) (2.2) 

The "unitarity method" H, |5| combines the unitarity cuts with the results of reduction to an expansion 
in master integrals Ii H 

A l-loop = ^ CiJ . (23) 

i 

The master integrals in d dimensions with massless propagators are scalar pentagons, scalar boxes, scalar 
triangles, and scalar bubbles. In the full d-dimensional formalism, there are no cut-free terms. 
The n-point scalar integral with massless propagators is 

~i(^) D/2 I 7^77 • (2-4) 



The coefficients Cj in (2.3) are, by construction, cut-free rational functions. In the unitarity method, 
we do not derive the coefficients of master integrals by performing any reduction. Rather, we take the 
coefficients as unknowns and proceed to constrain them by performing cuts on both sides of ( |2.3j ): 



Aj4 l-loop = ^- CiA/ ._ ( 25 ) 

i 

Any realization of the unitarity method must address the problem of isolating the individual coefficients 
q. The unitarity method succeeds because the cuts of master integrals are logarithms of unique functions 
of the kinematic invariants. 



In [13], it was shown how to obtain scalar box coefficients directly by cutting four propagators rather 
than two. Similarly explicit analytic formulas for the other coefficients have recently become available 

m n mm- 



Here, we refer to the formulas given in [27], after setting propagator masses to zero for simplicity. The 
generalization to the case of massive propagators has now been given in [46]. 

2.2 Input quantities 

Having present the general picture, in this subsection we can start with the following most general expres- 
sion for a unitarity cut integral: 

C= f d 4 -^c( / , 2 )ffiiL^l 5 ( + ) (p2)5 ( + ) ((p _^ )2) . (2 . 6) 



2 The delta functions here should properly be denoted by , indicating that they are restricted to the positive light cone. 
We shall drop the superscript for simplicity. 
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We work in the four-dimensional helicity scheme, so that all external momenta K{ are 4-dimensional and 
only the internal momentum p is (4 — 2e)-dimensional. We decompose the (4 — 2e)-dimensional loop 
momentum as 32] 



p = £ + ft, 



(2.7) 



where £ is 4-dimensional and jl is (— 2e)-dimensional. With the integrand in the form of ( |2.6| ), there is a 
prefactor c(/u 2 ) which depends on the external momenta as well as on ju 2 . In this discussion we shall be 
paying careful attention to all dependence on /x 2 . 



From this starting point, the coefficients of master integrals were listed in |27J]. Now, we would like to 
be able produce the complete 4-dimensional expression, by performing the integral over /i 2 by the recursion 
and reduction formulas of [18, 21 1. To get this complete answer, we need to consider the dependence of 
the prefactor c(fi 2 ) on fi 2 , along with the power of /i 2 in the coefficient formulas of p7f| . We consider this 
dependence in terms of the dimensionless parameter u, defined by 



u 



K 2 ' 



{21 



With this definition, the cut integral (^^) can then be rewritten as 



C 



r(-c) 



K 2 
4 



du u 



ifj 2 ) J (£ di) [£ d£] VT 



It 



(i\K\£] n+2 nil 



(2.9) 



The coefficients listed in [27], which are summarized below, are the results of the four-dimensional part 
of the integral ( |2.6D ; they are functions of u. The "four-dimensional cut-constructible" part of the amplitude 
could be obtained by setting u — > in each of these coefficients, inside the expansion of the amplitude in 
master integrals. The complete D- dimensional amplitude requires dealing with this u-dependence. Here it 
is enough to apply the polynomial reduction identities given in [18, 21|. These identities assume polynomial 
structure of the coefficients C(u), which is proven in the present paper, and which may also be deduced 
within other approaches |35[| . However, if we desire a result only through O(e ), it may be more efficient 



to use the dimensionally-shifted basis discussed in [49, 29 1. We shall return to this point in the following 
subsection. 

From the initial expression ([2j]), we extract all the necessary information, as follows. First, notice the 
triplet of integers 



(m, k, n = m — k) 



(2.10) 



which will play an important role. In particular, the value of n constrains the basis of master integrals 
P, P|- If n < —2, there are contributions only from boxes and pentagons. If n > — 1, contributions from 
triangles will kick in, and finally if n > 0, bubble contributions show up as well. This pattern is well known 
from traditional reduction techniques. 
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Second, we use the values of K, Pi, and Kj from the expression (^) to define the vectors Qj, Rj, and 
related important quantities, as follows: 3 



V^Kj-^K 
= K] — Kj ■ K 

Pi ■ K 

Qj(u) = — - u)qj + ajK, 

-2 



■WT^)Ki V^^^g- K 



Rj{u) = -(y/T=u) Pj + (3jK 

/ J 

K 2 



(VT^)Pj - (i - S^)ElJL k 



2.11) 
2.12) 
2.13) 
2.14) 

2.15) 
2.16) 

2.17) 
2.18) 



One important observation is that 

qj ■ K = pj ■ K = 0. (2.19) 
At this point, we wish to make a few more remarks. 

• The input quantities are given by K, Kj , Pj . From this we can define qj , ctj , pj , (3j and Qj (u) , Rj (u) . 
We make reference to the number of these vectors, encoded in the triple of integers (m, k, n). 

• To simplify notation when we set u = 0, we will write expressions such as Qj(u = 0), or just Qj. 

• The coefficients of the master integrals are polynomials in u. In this paper, we shall find that the 
maximum degrees of these polynomials are the following. Pentagon: 0. Box: [(n + 2)/2]. Triangle: 
[(n + l)/2]. Bubble: [n/2]. Here, [d] denotes the greatest integer less than or equal to d. 

For a renormalizable theory we have n < 2; thus we have the maximum degrees of 2 for boxes, 1 for 
triangles, and 1 for bubbles. These degrees are consistent with |L7], 22] and [29]. 



Knowing the maximum value of the degree of the polynomial in u, we can then calculate the coefficient 
of u s by the formula 

1 d s C(u) 



?! du s 



(2.20) 

u^0 



3 These definitions apply specifically to the case with massless propagators. Only a slight modification is necessary for 



massive propagators . 
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so 



1 d s C{u 

C(u) = > — — 

v ; ^ s\ du s 



s=0 



u . 



(2.21) 



The expression ( |2.20| ) is central in this paper. Since c s now has an expression where u does not 
appear (as indicated by the right-hand-side of expression ( |2.2C| )), it can be evaluated numerically. 4 

Summary of coefficients of 4-dimensional master integrals: 



For the box coefficient with momenta K, K r , K s , 

-2\2+n / Ttk+n 



C[Q r ,Q s ,K] 



Wj=l {Psr,l\Rj\Psr,2\ 



2 \(Psr,l\K\P sr ,2] n+2 Y$ =1 ^ (P sr ,l\Qt\Psr,2\ 

For the triangle coefficient with momenta K, K s , 
(K 2 ) 1+n 



+ {P* 



srA 



Psr,l\ 



C[Q S ,K] 



1 



1 



2 ( x /SI)"+ 1 (n + l)!(P S)1 P s , 2 



n+l 



d „ +1 / 0^(^,1-^,21^^1^,1-^,2) rD 

X~ — I — +{-Ps,l 



dr n+i ^J^^ (Ps)1 _ r P s , 2 |Q t Q s |P Sil - rP s , 2 
For the bubble coefficient with momentum iC, 

J n,n—q 



Ps,2} 



r=0 



= (^ 2 ) 1+ - E ( BgUoo +EE (fi££? 1} w - <^ ;2) ( S ; 



g=0 



where we have made the following definitions: 



1 



r=l a=g 



(2r? -Kr 1 n^f B (^(i£' + «7)W 



s=0 



dr n yn\[ V \rjK\n] n (t + 1)(K^ y ^ rg =1 (£|Q p (tf + sri)\i) 

l) b+l d b ( 1 {Pr,l-TP r MPr,l\ t+1 

lt+1 



T=0 



{Pr,i ~ TP r>2 \Q r r]\Pr,i ~ TP r>2 ) b l\]=i (fi.i ~ rP^R^K + sr?)^ - rP r , 2 ) 
(P r ,i - rP r!2 |r/K|P r ,i - rP r / +1 ]lj=i, p ^ < P r,i - rP r ^\Q p {K + s V )\P r>1 - rP r , 



T = 



(-1) 



6+1 



d b ( 1 <Pr,2-rP r ,i|7/|P r , 2 ] 



t+1 



6!^A; b+1 (Pr,l Pr,2) bdTb \( t + 1 ) (P r ,2-rP rjl \K\P r>2 } t+1 
(P r ,2 ~ rP r>1 \Q r ri\Pr,2 - rP r ,i) b U]=i ( p r,2 ~ TP r ,l\Rj(K + srj)\P r , 2 - rP r ,i) 



(Pr,2 - rP r ,i|77^|P rj2 - rP rjl ) n+1 nj=i, p ^ ( p r,2 - rP r ,i|Q p (K + Sr))\P r>2 - tP, 



T=0 



Note that the prefactor c(^ 2 ) has not been included in these formulas for coefficients. 



4 See [M for another approach that is possibly more efficient. 
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2.3 Some important constructions from input quantities 

Given two momenta S, R, we construct two null momenta. If R 2 = and S 2 = 0, R,S are themselves the 
two null momenta. If at least one of them is not null, for example R 2 ^ 0, then we construct two null 
momenta as follows. 



-2S ■ R + JA(S, R) -25 • R - JA(S, R) 

P(S,R);i = S + X * R > X ± = 2fi 2 - » *2 = ^ 2 y \ (2.22) 

where 

A{S,R) = (2R- S) 2 -4R 2 S 2 . (2.23) 
Then, the following quantities necessarily vanish. 

= <P(s, Ji) ;i|S r |i'(5,fl). 2 ] = (P(S,R);2\ S \ P (S,R);l] = (P(S,RyM R \ P (S,R);2} = (P(S,R);2\R\P(S,R);l] • (2-24) 

We shall use the following identity: 

(Pl\V\P 2 \ {P 2 \W\P{\ = tr (^—^- Pi y fa W 



= i((2Pi • V){2P 2 ■ W) + (2Py ■ W)(2P 2 ■ V) - (2P X ■ P 2 )(2V ■ W) - Ue^ ap P? V V P%W>). (2.25) 

Any four-dimensional momentum K can be expanded in a basis of four other independent momenta 
Ki,Kj,K m ,K n by 

K = a m K m + a n K n + ciiKi + ajKj, (2.26) 

where the coefficients are given by 

t(Ki,K, Kj,K n ) e(Ki,K m ,Kj,K) 

dm 



m e(Ki,K m ,Kj,K n y n t{K h K m ,K h K n y 

= e(K, K m , Kj,K n ) = e(K u K m ,K,K n ) 

ai e(K{, K m , Kj,K n ) ' aj e(K{, K m , Kj,K n ) ' 1 " ) 



with 



e{K u K 2 , K 3 , Ki) = e, lvpi K{K v 2 K p z K\. (2.28) 
At times, we will write these coefficients in the form a ^> K i> K ™' Kn > K ) to emphasize the related quantities. 
2.4 The dimensionally shifted basis 

Since, as we shall demonstrate, our coefficients are polynomials in u, we can translate this information into 
the dimensionally shifted basis [49]. More explicitly, if we define 

f d D P P ai p a -m 

I D \P a \..P a ™} = -iUir) D/2 / y-4n 1 . ( 2 - 2 9) 

J (2vr) B (pa _ m 2)...((p -YTiZlh) 2 -m 2 ) V ' 
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then we have 



lt 2e V) k ) = ^z^I ( rt 2e+2k) [l]- (2-30) 

There are two merits of using this basis. First, we can throw away all 0(e) contributions to make 
the calculation easier. Second, we improve efficiency. To use the recursion and reduction relations, we 
first calculated all the contributions by reduction to boxes, triangles and bubbles, and then added them 
up. With the dimensionally shifted basis, this process of reduction/summation can be done in one step, 
simplifying calculations. The usefulness of this dimensionally shifted basis has been discussed in [22, 29]. 
Here, for reference, we discuss this evaluation in Appendix |A[ 



3. Triangle coefficients 

Now that we have the necessary background information, it is simplest to start with the coefficients of 
triangles. Some features of this discussion will apply to bubbles as well. 



3.1 Simplifying the formula 

We write the formula for triangle coefficients from 
u-dependence. 



C[Q s (u),K] 



(K 



2\l+ra 



in the notation of the previous section, emphasizing 



(3.1) 



v n+l 



2(VA(Q s («),if)(n))«+i ( n + 1)! (P {Qs{u)>K );i(u) P(q,(„),jc); 3 (« 

d n+1 ( \\)ti (P{Q s (u),K);l( U ) ~ TP (Qs(u),K)A U )\ R j( U )Qs( U )\ P (Qs(u),Ky,l( U ) ~ tP (Q s (u),K);2( U )) 

X "~~ ^ I T 

dT " Vn t =Ws ( P (Q s (u),K);l( U ) ~ TP (Qs(u),K);2( U )\Qd U )Qs( u )\ P (Q s (u),K);l(u) ~ tP(q s{u)]K) . 2 (u)) 
+{ P (Q s (u),K);l( U ) P (Q s (u),K);2{u)}) 

/ T — 

where the P(q s ( u ),k);1,2( u )> as depicted in the indices, are constructed in terms of Q s (u),K, as defined 
in ( |2.22|) , and depend on u. In principle we can put (3.1) into ( |2,20| ) to take derivatives. However, the 



u-dependence everywhere might be an obstacle to taking stable derivatives in terms of u in ( 2.20 ). In this 
subsection, we recast this u-dependence in a simpler form. 

Using the definition of Q s from Q2.15 ), and the property ( [2.19 ), we find from ( |2.22 ), ( 2.23j ) that 



A(Q s (u),K)(u) 



(l-u)(-4q 2 s K 2 ), 



-2a s K 2 ±^A(Q s (u),K) 
2K 2 
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When we take the square root of A(Q s (u), K), there is a sign ambiguity. It can be shown that the choice 
of sign does not affect the final result. To be explicit, we choose the minus sign here, i.e., 



y/A(Q 8 (u),K)(u) 
xi t2 (u) 



VT^l^-4q 2 K 2 



K 2 



-Vl -« 

= -(Vl -«)yi,2 - ol s , 

where we have defined new scalar quantities, y\ 2, as follows: 



a. 



J{K S -K) 2 -K 2 K 2 
2/1,2 = ± — = ±- 



if 2 



K 2 



With these results, we can see that 



P{Q s {u),K);i{u) = - (Vl - u)q s + a- s # + Xi(u)K = -y/Y^u(q s + = -(Vl - u)P { 



(3.2) 
(3.3) 
(3.4) 

(3.5) 
(3.6) 



The u-dependence has been factored out; here the null momentum P( qs: K);i does not depend on u, since it 
is constructed from q s , K-as indicated in the subscript indices. 

Substituting ( |3.2| ) and (^) back into ([O]), we find that the factor Vl — u h as cancelled out. Thus 
we have 

(K 2 ) 1+n 1 



C[Q S ,K] 



jn+l 



2(_Vr^)«+i(v / A(CT) n+1 (n + l)l(P M;1 P W! ) 

IIj=l ( P (g a ,A:);l - ^(g^^^l^jNQsNI^,^);! - tP^ k) . 2 ) 

dT n + l [ U k t=1Ms (P {qs ,K);l - rP {qs , K) . 2 \Qt(u)Q s (u)\P iqsjK y A - tP^ K) , 2 ) 

To simplify further, apply the identity (t\QQ\t) = to derive 



n+l 



+ {P(q s ,K);l P (qs,K);2} 



{£\Qt(u)Q s (u)\£} = (£\(Qt(u 



Q s («))Q a («)|^ 



-VT^A|fe--g s )Q,( 

\ a s 



u)\l 



Rj{u)Q 8 {u)\e) = U\{Rj{u) - ^-Qs{u))Q s \£ 



Pi 



If we define two more vectors qt , Pj by 



at 



qt = (qt q s ), 



-VT^(£\(p j -^q s )Qs(u)\£ 

Ota 



Pi 



Pj = (Pj -Is), 



(3.7) 
(3.8) 

(3.9) 



then we use the identities (|3.7|), ( |3.8[ ) to conclude that 
(K 2 ) 1+n 1 



C[Q 8 ,K] 



"(q a ,K);l P (q s ,K);2 



n+l 



2(VA(g s ,^))"+ 1 (n + l)!(P (( 
dr n+1 ynLi.t^ < P (g s ,^);i ~ rP (g s ,K);2l5tQs(u)|P ( g ai ^ ); i - tP^ k) . 2 ) 



+ {-Pfe.E-);! ^ P (q a ,K);2} 



(3.10) 
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Compared to ( |3.1| ), the u-dependence in ( |3.10| ) is much simpler; all u-dependence here comes only from 
Q s (u). Thus, ( |3. 1Q|) is well suited for use in ( p. 20 ). 



How to use the formula ( |3.10| ): The degree of this polynomial in u will be seen to be [(n + l)/2]. 
Thus we can get the corresponding coefficients by taking derivatives in u first (from to [(n + 1)/2], to get 
coefficients from each term in the polynomial), and then setting u = 0. 

For example when n = — 1, we can set u = directly and get 



C[Q S ,K] 

ne{o,-i} 



(K 



2\ra+l 



1 



n+l 



d n+1 ( Yijti {P{ qs ,K);i ~ rP {qM \pjQ a \P {qM - rP (qs , K) , 2 ) 
dT n + i \^ ( P(?sjK);1 _ rP (fc>J0 . a |ftQ.|P (fcfJOsl - rP { 



+ { P (q a ,K);l ^ P (q s ,K);2} 



(3.11) 



which is suitable for numerical evaluation. For n = 1, 2 the result will take the form of a linear polynomial, 
cq + c\u. To get ci we take one derivative, using 



dQ s (u) 



du 



u=0 



2 



(3.12) 



In Appendix [b], we have explicit expressions, free of derivatives, for triangle coefficients when n < 2. 
The formula (3.10) contains u in both numerator and denominator, so it is not so obvious that the 
total result is simply a polynomial in u. The proof of this property is given in the next subsection. 

3.2 Proof that the triangle coefficient is a polynomial in u 

We start by considering two quantities that arise in our expressions, in the course of taking derivatives: 



El = (P( q s,K);2\PjQs(u)\P( qs ,K);l) + ( P (q s ,K );1 \PjQs(u) |P(g S) iC);2> 
^2 = {P{ qs ,K);l\PjQs{u)\P {qs ,K);l) ( P (q s ,K );2 \PjQs{u) |P(g 3 ,X);2> 



By writing Q s (u) as a linear combination of the P( qs; K);ii 



Qs(u) 



+ 



P 



2 2yi 

and recalling that g s • K = pj ■ K = 0, we find that 

a.K 2 



Vl-u Og_ 
2 2yi 



P 



(9s,^);2' 



(3.13) 
(3.14) 



(3.15) 



All n-dependence has dropped out of this expression 



(2pj ■ q s ) (P ((?Sij k) ; i P( q3y K);2) ■ 



(3.16) 
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For E 2 , similar manipulations show that 

E2 = (P( qs ,K);l P( qs> K);2) 2 ~ (Qs • Pjf) + 1 - «) , (3.17) 

which is a polynomial in u. 

Now we prove that the full expression ( |3.10| ) for the triangle coefficient is a polynomial in u. Throughout 
this proof, let us abbreviate Q s by Q and P( qst K);i by Pi- 

The triangle coefficient is given in terms of derivatives with respect to t on an expression where the 
r-dependence appears in the factors (Pi — TP2\PjQ{u)\Pi — TP2) (in numerator or denominator). After 
taking the derivatives, we set r = 0. In this process we will produce the following three combinations: 
(Pi\pjQ(u)\Pi); E\\ E 2 . It is easy to see how the first two combinations arise. The third combination, E2, 
appears, for example, in 

(P (Pi - rP 2 \p j Q(u)\P 1 - tP 2 ) _ -2 (P^pjQju)^} (P 2 \p j Q(u)\P 2 ) 
dr 2 (P^pjQiu)^} 

Consider the r-dependent factors in the denominator. With each derivative, we effectively add 
one overall factor of (Pi — tP2\XQ(u)\P\ —TP2) in the denominator and place one new factor, either 
{{P 1 -tP 2 \XQ{u)\P 2 ) + (P 2 \XQ(u)\P 1 -tP 2 )) or (-2 (P 2 \XQ(u)\P 2 ) (Pi|XQ(u)|Pi» in the numerator. 
After taking n + 1 derivatives, there are (n + 1) additional factors (Pi — tP2\XQ{u)\P\ — TP2) in the de- 
nominator. Thus, we have exactly k + n factors of (Pi\XQ(u)\Pi) in both numerator and denominator 
(after taking the r — ► limit). The u-dependence is exactly cancelled in this part, since we have 

(Pi[XQ^)|Pi) _ (P|A|P 2 ] 

(PilX'QHlP) (Pi|A'|P 2 ]- 1 ' ; 

The remaining u-dependence comes only through the factor E 2 in the numerator. By our previous calcula- 
tions (HH) and ( pTD , we see that indeed our final expression is a polynomial in u. Since every sequence 



of two derivatives in ( 3.1C| ) will produce one E2 factor, the degree of the polynomial is [(n + l)/2]. 



4. Bubble coefficients 

Our proof that bubble coefficients are polynomials in u is more complicated than the one for triangle 
coefficients, so many of the details have been relegated to Appendix p. Here we present the simplification 
of the li-dependence in the formula for bubble coefficients. The idea of the proof is consider the bubble 
coefficients are polynomials in \f\ — u, and show that the odd powers vanish. 

4.1 Simplification 

The coefficient of the bubble integral with momentum K is given by |27]] 



C [K] = ( K r +n E (fi£l>) + E E w - w 

q=0 \ r=l a=q 



(4.1) 

s=0 
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where 



B 



(«) = 



<f / 1 (2r]-K) t + 1 

y n \[ V \7jK\ V } n (t + l)(^ 2 ) 4 +i \t ^"+1 n j =i (e\Q p ( u )(K + sr,)\£) 



♦ |Jf-T^l7(] 



(4.2) 



T=0 



6! v /A(Q I .(w),ii') b+1 (P(Q r („),K) ; i(u) P(Q r («),ii:);2(w)) 

(P(Q r (u),K);l(u) ~ T--P(Q r (u),if);2(u)|?7|P(Q r ( u ),ii:) ; l(u)] 



t+1 



( + 1 



dr" \(t + l) (i^ gr(u)iJr) . 1 (t»)-rP ( g r(u)jJOi a(«)|ii'|i , (Q r («),JOii(«)] 

- ^P(Q r ( M ),iQ;2(")IQr(M)77|P(Q r ( M ),K);l(") ~ tP(q,, (m) ^) ;2 (m)) 
(P(QAu),K);l(u) ~ TP(Q r(u))J i ir );2(«)l»7-K'|-P(Q r (u),K);l(u) - TP(Q r («),K);2("))" 

II"^l fc (P(Q r («),K);l(") " TP(Q r ( M ),K);2(M)|^j(M)(^ + g??)|P(Q,(u),g);l(M) ~ T%( u ), 

X T 

E[p=l,p#r (P(Q r .(«),K);l(w) - rP(Q r(u)) j ir );2(w)|Qp(«)(-K' + S??)|P(Q r («),if);l(w) - TP (Qr{u y K) . 2 (u)) / 

W v /A(Q r (u),^) (P(Q r („),K);lW P(Q r («),K);2(w)) 

rf b / 1 ( f (Q r (u),g);2(«) ~ rP {Qr(u)7K)ll (u)\r]\P (Qr(ulK) . 2 (u)] t+1 
dr» \(t + l) (P {Qr{uhKh2 (u) -rP (Qr(u) ^ );1 ( W )|^|P (Qr(u)iX);2 ( M )]* +1 

(P(Q r (u).K);2(M) ~ TP(Q r .( M ),K);l(")IQr(M)?7|P(Q T .(n),A:);2(") ~ ^%(n),g);lM)'' 

x ^ i ^ 

(P(Q,»,K);2<X) ~ 7-P(Q r («),Jif);l(«)|»7^|P(Q r (u),X-);2(«) - 1"P(Q r («),*T);l («)) 

(PtQ.W.^^M ~ TP^^^gHlPj^)^ + SV)\P(Q r (u),K);2( U ) ~ ^P(Q,(«),if);l( U )) 



(4.3) 



r=0 



• (4-4) 



T = 



n p =l,p#r (P(Q r (u),K);2(u) - TP(Q r ( u ),K);l(u)\Qp(u)(K + sr])\P {Qr{u) , K) . 2 {u) - T P {Qr ( U ),K);1 («)) , 

In this expression, since P(Q r ( u ) ; _ft-) : i('u) and A(Q r (u), K) depend on u, we want to simplify the above 
expressions as we did in the triangle case. Using (|3.2| ) and ( |3.6D , we see that we can pull out some factors 

5 



of VT 



u, giving 



B 



(r;6;l) 



1 



n,t 



S, U 



b\{VT^) b+ WHlr,K) h+l (P (qr , K) .l P iqr ,K);2) b 
1 ( P (qr,K)-A ~ rP {qr:K) . 2 \v\P(gr,Ky,l] 



1 1+1 



dr» \(t + l) ( P _ rP (grii , );2 |K|P (grii , );i ; 



t+l 



- TP( qr ,K);2\riK\P {qr:K y A - TPfajc). 



, n+1 



I\]=i (P( qr ,K)-,i ~ rP^y^Rjiu^K + srj)\P {qr>K) . x - TP { 



(qr,K);2, 



Up=i,p^r ( P (gr,K);i ~ rP {qrtK) . 2 \Q p (u){K + sri)\P {qriK) . x - TP [qr>K) . 2 ) 



(4.5) 



T=0 



5 One way to see it is to choose \P(Q T ( u ),K);i) 
cancels out immediately. 



\P(q r ,K);i) and \P(Q r ( u ),K);i] 



— u\P(q ri K);i]j the factor —yl — u 
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Similarly, for the function B^f (s,u), 



Bnt(s,u) = — — — , 

b\(y/T^H)^y/A(q r ,K) + (P {gr ,K);l P( qr ,K);2) 
d b ( 1 {P( qr ,K);2-rP {q ^ K )MP(Qr,K);2] t+l 



1 



t+l 

b 



^ V t + V (P (qr ,K)-2 ~ TP {9rtm \K\P (9rtK);1 

( P (qr,K);2 ~ ^P(q r) K);l\Qr{u)ri\P {q ^ K) .2 ~ TPj^K)^ 

X n~\-l~ 

(P(q r ,K);2 ~ TPfaMtllvKlPfcMp ~ tP^K)^ 

II"=1 (P(q r ,K);2 ~ ^P(q r ,K)A\Rj{u){K + Srj)\P {q ^ K) . 2 - TP^ rjK);1 ) 

Up=l,p^r ( P (qr,K);2 ~ rP(q r , K );l\Qp(u)(K + ST])\P {q ^ K) . 2 - TP {q ^ K) . l/ 



(4.6) 



T=0 



So, the bubble coefficient is now given by ( |4.1| ), ( |4.2|) , ( |4.5[ ) and ( |4.6| ). Unlike the triangle, since we 
now have factors of the form (£\Q(u)(K + sr/)\£) instead of {£\Q(u)K\£) , we cannot pull out further factors 
of \/l — u. In these expressions, the only u-dependence is coming from Rj(u) and Qt(u), which is much 
simpler than the original formula. 

The best way to use these formulas is similar to the triangle case. We can see from the formulas ( [4.2[ ), 



(4.5), fl4.6|) that the degree of the polynomial in u is [n/2]. Thus we can get the corresponding coefficients 



by taking derivatives with respect to u and then setting u = 0, as in ( |2.20|) . 
5. The box and pentagon coefficients 

Now we consider the box coefficients. This is the most complicated part, although the formula is the 
simplest! There are two reasons for the complexity. First, the box coefficients contain not only true box 
coefficients, but also pentagon contributions, indicated by a linear factor (au + b) in the denominator. 
We should be able to separate the box part from the pentagon part. The second reason is that the null 
momenta Piq^^.q^u))-^ depend on u in a very nontrivial way (as Qj{u) + x a Qi(u)), unlike the cases of 
triangles and bubbles. 

Given the vectors Qi,Qj,K we can construct a vector (^■" q 3' K ' orthogonal to all three: 



(qot' q - K) = ± e ^q^ (5.1) 



K 2 

We shall make use of our collection of orthogonality relations: 



^ vpi K\K^. (5.2) 



qi -K = q r K = q -K = (5.3) 

qo ■ m = qo ■ qj = qo ■ K = o (5.4) 
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Observe the fixed ordering of qi,qj in the definition (5.2) of 



(qi,qj,K) 



Exchanging them leads to a minus 



sign difference. The ordering is connected with our definition Pji = Qj + xQ{. 

Because new features arise in different cases, we have divided our discussion into the three cases k = 2, 
k = 3, and k > 4. 

In the process of simplifying our expression, we also demonstrate that the coefficients are truly poly- 
nomials in u. Thus this is a constructive proof. 

5.1 The case k = 2 

In this case, there is exactly one box, whose coefficient is given by 



(K 



2\2+n 



( IIS ( P (Qj («),Qi(«));l ( n ) I R s ( n ) I P {Q 3 (u),Qi(u));2 (u) 



\ 



\ \ i> Wi(«),Qi(«));i( u )l- K 'l- p (Qi(«),<3i(«))j2(«) 



n+2 



+ { P (Qj(u),Q t (u));l( U ) ^ ^(QiW.QiH)^^)} 



/ 



Later, we will reduce the other cases to this expression as well. Now we carry out a detailed calculation to 
show that the u-dependence is polynomial, and we find an expression where the u-dependence is easier to 
see. 

In the numerator, the factors R s (u) have ii-dependence in the form R s (u) = —{y/l — u)p s + S K (here 
we do not assume any particular form of /3 S ), with p s ■ K = 0. Thus the vector p s can be expanded in a 
basis of the three-dimensional vectorspace orthogonal to K, as follows. 



Ps = % 

The coefficients in this expansion are: 



(qi,qj,K;p s ) (quqj,K) (qi,qj,K;p a ) {qi,q.j,K;p s ) 



{q uqj ,K; Ps ) _ (Ps • Qihj ~ (Ps ■ fcX?* ' Qj) 



qh] - (ft • ijf 



{9i,qjJC;p.) = (Ps ; gj)g, 2 ~ (Ps ; gi)(gt • Qj) 

3 qh] - (« • 

( qi , qj> K;p s ) _ (p s ■ qQ l,q3,K) ) _ e{p s ,q i ,q j ,K) _ e(P s , Kj, Kj,K) 



Using this expansion, we can write 

R s ( u ) = a^> K]Pa \-(VT^)q ( t qj ' K) ) + af i,qhK ' pa) Qiiu) + af^^ Pa) Q 3 {u) + $»«> K ^K, (5.9) 
where we have defined 

Using the result fl|gg ) that 



(qi,qj,K), 



(qi,qj,K)\2 



(5.5) 

(5.6) 
(5.7) 

(5.8) 



^(Q i (u),Q i (u));l(«)l<9i/i( n )l jP (Qi(«) ! Qi(«));2(' lt ) 



0, 



(5.11) 
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we have 



jP (Q J H,Q l («));i( u )l jR s(u)|P(Q j ( u ) i Q i ( u )) ;2 (n 

- — ( y/l — ,,\J qiS ^ K ' P ^ 



u a r 



P{Qj{u),Q i (u))-,l{u)\q^ l ' q3 ' K \P(Q j ( u ),Q i {u))-2{ u ) 



Thus, we can write 



Us=i (■ p «j(u),0 i («));l( u )l- R »( u )l- p (0 J («),0i(«));2(«) 



- P (Qi("),Q l («));i( M )l i;C l P (Q J («),Q l («));2( M ) 



n+2 



(5.12) 



(5.13) 



= £3 

/i=0 

Here we have defined 



n+2 (g<iq . |J0 (-V^^ fe (^(Q>)^^^ 

hi -\h 
h " < P(Q,(u),Q,(u));l{u)\K\P {Qj{u)yQi{u)) . 2 (u 



a 



(qi,qj,K) 
h 



e n- 

SC{l,2,...,n+2} SGS 

\S\=h 



(qi,qj,K;p s ) TT R {qi,qj,K;p a ) 



(5.14) 



ses c 



where 5 C denotes the complement of S: S c = {1, 2, . . . , n + 2}\5. 

Now we can show that the box coefficients are indeed polynomials in u, in this case where k = 2. By 
the above expansion, the coefficients are given by sum of the following typical terms. (To simplify the 
formulas in this proof, we shall now write Pji >a (u) in place of P(Q j ( u ),Q i {u)):a{ u )-) 



(P jit i(u)\K\P jifi (u)]' 



+ 



(P ji , 2 (u)\K\P ji>1 (u)] h 



(5.15) 



( 



(P i i 1 l(«)|^|Pi^(«)] & (Pii,2(«)|^|Pii,l(u)]' 1 



(-VT^) h (2t-(l - u) y /^{Q i {u),Q j {u))qlK 2 ) h + (-2i(l - «) VA(Qi(u), Q»)ggtf 2 )' 



(5.16) 



if2(l- u )2[( 2 g i .^) 2 -4g 2 ^] 

We have used the following results (making repeated use of the identity ( 2.25| )) 

K 2 



(P jiA (u) \K\P jij2 («)] (P iii2 («) | K | P ji)X («)] 



Qi{uf 



{l-uY[{2 qi .q 3 Y-Aqiqj\ 



{PjiM\to\Pji,M\ (P 3 i,2{u)\K\P Jhl {u)} = 2i(l-u)^ ^f:® j{u)) q 2 K 2 



<P ji ,i(«)|lif|Pii ) 2(«)](^i,2HI?o|i'ii,i(«)] = -2»(1 ~ n)^ ^^ 2 9j(M)) g0 2 ^ 2 
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and 



A(Qi(u), Qj(u)) = (1 - u) {(1 - u)[(2g* • </,) 2 - 4g?g|] + 4K 2 [a 4ai (2 % • ^) - a 2 ,? 2 - a 2 g 2 ]} . (5.17) 



Equation ( 5.16 ) is our most important result in this subsection. There we can see that when h is odd, we 
get zero. When h is even, the factor (1 — u) 2h in the denominator is cancelled by corresponding factors in 
the numerator. (Notice the overall factor of (1 — u) within A(Qi(u), Qj(u)).) Thus we see that indeed the 
expression ( |5.15D is a polynomial in u, specifically, 

[{2qi . qj) 2_ 4qh , ]h tor n even, ^ lg) 

for h odd. 

It is clear that the maximum degree of the polynomial is [^p]- 

5.1.1 A simpler expression for k = 2 



The aim of this sub-subsection is to find another expression with the same value as ( 5.18| ), but with more 
transparent u-dependence. We have just seen that all li-dependence has cancelled out, except for the 
second factor of (1 — u) in the first term of Ay(u), as it appears in ( |5.17| ). Since Q5.18| ) is now an expression 
in terms of scalar quantities, we can consider the effect of setting u = at the beginning of the calculation. 
Recovering the expression ( 5.18| ) then requires the following single modification: 



[a iaj (2 qi .qj )-q?qf - qffi] 



AK 2 [a i a j (2q i -q : j)-a 2 q 2 -a 2 q 2 ] 



1 + 



(2q l -q J ) 2 -4q 2 q 2 



Then, 



(-^T=^)h (P jiA ( u )\q \P jit2 (u)] h {-^T^u) h (Pji, 2 (u)\qo\Pji,i(u)} 



+ 



(P jiil (u)\K\P jit2 (u)] h (P ji , 2 (u)\K\P j i >1 (u)] h 

_ {P 3l ,i{u = 0)\a^(u)(- qo )\P 3l!2 (u = 0)} h + (P ji>2 (u = 0)|«fa^)(n)(-g )|P^ 1 (^ = 0)]" 
(Pji^u = 0)\K\P jij2 (u = 0)} h ' (P jh2 (u = 0)\K\P jiA (u = 0)} h 

Now, all of the u-dependence is concentrated within the a^ guqj \u). Going back to ( |5.13| ), we can perform 
a similar operation: 

Uti (^Au)\R s (u)\P^ 2 ( U )] 

^+2 1" i^yW Pji,2{ u )S 



(P jiA (u)\K\P jii2 (u)] 

n:i 2 (pjuiu = o)\R s ( U )\p jit2 ( U = o) 



(P jitl (u = 0)\K\P jit2 (i 



in+2 



+ {P jiA (u = 0) <-> P jii2 (u = 0)}. (5.21) 
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where we have defined 

R s (u) = a { ^' K;Ps) a^\u)(- qo ) + pi q ^' K]Pa) K. 
With the definition ( pM0|) , we can rewrite R s (u) as 



(5.22) 



R s = a ^' K ^\a^(u) - l)(-g ) - a^'^Q^u = 0) - af"^ K ' Ps) Q 3 (u = 0) + = 0). 



Now we make use of the properties 



(P jiA (u = 0)\Qi(u = 0)\P jit2 (u 



(P Ji: i(u = 0)\Qj{u = 0)\Pjw(u = 0)] = 0, 



and the fact that R s (u) is defined such that the equation ( 5.2 1] ) is satisfied, we can drop the terms with 
Qi/j{u = 0), so that 



Rs(u) 



ps q ^ K \^ q . ){u) _ i){ _ q ^ K)) + Rs{u = Q) _ 



(qi,qj,K), 2 



(5.23) 



Equations <^2^) and (|3|) are our final simplest result. All n-dependence has been packaged inside 
(u) -1), which is zero when u = 0. Also, it has now become clear that the degree of the polynomial 
in u is [(n + 2)/2]. 

Summary: The box coefficient for k = 2 is given by 



C[Ki,Kj] k=2 



(K 



2\2+n 



' IIS ,00;! l-R. ( u ) I ,Qi);2 



V 



P (Q J -,Qi);il- K 1- P (Qi,Qi);2 



n+2 



with 



3. (it) 



(qi,qj,K) 
% (Aquqj) 



{qi,qj,K) y 



(a^>(u) - l)(-qi qi ' qj ' K) ) + R s (u = 0). 



(5.25) 



Let us emphasize again that P(Q j ^Q l )- a is constructed from Qj(u = 0) + x a Qi(u = 0), so it is independent 
of u. 



5.2 The case k = 3 

In this case, we will see the pentagon show up, and we shall learn how to separate boxes from pentagons. 
Again, we shall abbreviate the notation of the vector P(Q j ,Q i )- a by Pji,a- We evaluate an expression of the 
form: 



n:=! (p^)\R s (u)\p jh2 (u) 



(P jiA (u)\K\P ji)2 (u)] n+2 (Pji,i(u)\Qt(u)\Pjw(u)] 



+ {Pji,l(u) <-> Pji,2(u)}. 



(5.26) 
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Again, we would like to expand R s (u), or equivalently p s , in a suitable basis of the vectorspace orthogonal 
to K. In this case, we do 
qi,qj,qt available already. 



to K. In this case, we do not need to construct the vector q^ z ' q: " K \ because we now have three vectors 



p s = a [ qi ' qjmPs) q t + af- qM) qi + a {qi ' qjmPe) qj. (5.27) 



Thus, 



{P 3hl (u)\R s {u)\P jh2 {u)\ = a ^ qt;Ps) (P ji;1 (u)\Q t (u)\P jii2 (u)] + /^*^ (P Jhl (u)\K\P Jh2 {u)] , (5.28) 
where we have defined 

= (Aj _ J- a iq ^' qt ' Pa) a h ). (5.29) 

h=i,j,k 

Now we can expand the rational function within the coefficient formula: 
<P ji , 1 (u)|if|P ji , 2 (u)] n+2 (Pji^lQtiu^P^iu)} 

= (Pji,l( u )\Qt(u)\Pji,2(u)] h 3Q 

where we have defined 

i(<n,qj,qt) _ tt (qi,qj,qt;ps) tt o(qi,qj,qt;p s ) 



SC{l,2,...,n+3} sgS sGS c 
|S|=h 

Now, break the sum ( |5.30| ) into two parts, by separating the term with h = from the rest. In every 
term with h > 0, the factor {Pji^(u)\Qt{u)\Pji :2 {u)\ from the denominator is cancelled by the numerator. 
What remains is a term in the form we considered in the previous subsection, the case k = 2. That part 
contributes only to boxes. We shall return to that part in a moment, to find the exact box contribution. 
Notice that we can observe at this point, from the comparison to the k = 2 case, that the degree of the 
polynomial is again [(n + 2)/2]. 

The term in ( |5.30|) with h = is 

i(qi,qj,qt;ps) 

{P Jhl (u)\K\P ji<2 (u)} 



(PiiiMiQtMiPjiMY 



The cut of a pentagon integral has been analyzed in [21]. We clarify its u-dependent behavior in Appendix 
It is directly related to the sum of three cut-boxes. The part of the cut-pentagon that is related to the 
cut-box C[Qi,Qj] is 

{P^ x {u)\K\P jii2 {u)\ \ 

+ {P jitl (u) <-> Pji, 2 (u)\ . (5.31) 



2K 2 \(P jiA (u)\Q t (u)\P ji72 ( U )) 
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Thus, we see that this h = term is exactly a pentagon contribution. The coefficient of the pentagon 
integral must be 



C[Qi,Qj,Q t ] = (K 



n+3 

2\3+n TT «fe'9j.?t;Ps) 



(5.32) 



which is entirely independent of u. 

Now let us return to the box coefficients, using our result (5.24), along with the definition (5.23) 
applied here to the vector Qt(u) instead of R s (u). 



(<n,qj,Qt) 



h=l 
n+3 

h=l 

n+3 

h=0 
n+3 



(Pji^lKlPjiM] 



h-l 



P ji , 1 (u = 0)\Q t (u)\P ji:2 (u = 0) 



h-l 



(P ii ,i(« = 0)| J K'|P ii ,2(« = 0)] 



h-l 



P ji>1 (u = 0)\Qt(u)\P ji , 2 (u = 0) 



(P jitl (u = 0)\K \P ji>2 (u = O)]*- 1 {Pji^u = 0)\Q t (u)\P jit2 (u = 0) 
( qi , qj ,m P s) (P^i(n = 0)|^|P Jt , 2 (n = 0)] 



=i 



P jitl (u = 0)\Q t (u)\P ji , 2 {u = 0) 



n:i 3 ( Pf<,i(« = o)|p s (u)|p ji!2 ( W = o) 



(Pji.iCu = 0)\K\P ji>2 {u = 0)] n+2 (Pji^u = 0)\Qt(u)\Pji,2(u = 0) 



{qi,qj,qf,Ps) ( P ji,l( u = 0)\K\Pjj, 2 {u 



n+3 

7=i {PjiAu = o)\Qt{u)\Pji, 2 {u = o) 

In the last line of this equation, we have defined 



(5.33) 



l («,?,,« a ».) o («*,^») (a ( B> 5 J 0( u ) _ i)(_ go) + # s ( n = o) - J2 



(qi,qj,qt;ps) 



QJu = 0). 



1=h3 



The cumbersome double-tilde notation is a temporary inconvenience, to avoid confusion with R s (u), which 
was defined in (5.22) for the case k = 2. In fact, we shall discover later that the two quantities are identical. 
Using the property \Pji t i(u = 0)\Qi/j{u = 0)\Pji i2 (u = 0)] = 0, we can redefine this vector as follows: 



R s (u) = 4 q ^' qt ^ ] a^' K -' qt \a^\^ - l)(-q ) + R s (u = 0). 
Finally, the box coefficient is given by 



(5.34) 
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C[Qi, Qj]k=3 



(K 



2\2+n 



Il?=l ( P {Q 3 ,Qi);l\ R 3i u )\ P {Q 3 ,Q&,2 



n+2 



n+3 



P {Qj,Qi);x\ K \ P {Qj,Qi);2 



«=i (- P (Qj,Oi);il < 3*( u )l- P (Oj,0<);2 



P (Qj,Qi);l\Qt( u )\ P (.Q3,Qi);2 



(5.35) 



As an alternative to the expansion ( |5,27j ) of p s in the basis qi,qj,qt, we consider another expansion: 



E 



(K,Ki,Kj,K t ;P s ) jy- 



(5.36) 



By projecting equation ( 5.36 ) onto the vectorspace orthogonal to K, and comparing with ( 5.27 ), we see 
that 



{K,Ki,Kj,Kt;Ps ) _ (qi,qj,qt\Ps) _ ■ ■ . 

u w — u w , w — t, J, L. 



(5.37) 



The advantage of the expansion (5.36) in four vectors is that we can solve for the coefficients explicitly, 



as in (2.26), and find 



(q l ,q J ,qf,Ps) _ (K,Ki,Kj,K t ;P s ) = e(Ki,Kj,K,P s 

h ~ at ~ e{Ki,Kj,K,K t 



(5.38) 



Thus, using (|5.29[) , we have 

Mam**) _ K i< P *> K i' K > K ^ + K H K ^ P ^ K ' K t) + K 2 e(K i: Kj, P„K t ) + K?e(Ki, Kj, K, P.) 



K 2 e(Ki, Kj,K, K t 



With the formulas ( |5.3&j ) and ( |5.8| ) for the expansion coefficients written in terms of input vectors, we can 
simplify the following coefficient in the definition ( 5.33; ) of R s (u): 



(q l ,q J! qt\Ps) (q t ,q 3 ,K m ) __ e(Kj, Kj , K, P s ) e{K t , Kj , Kj , K) _ e{P s ,Ki,Kj,K) _ ( qi ,qj,K;p a ) 



We conclude that 



e(K t ,Ki,Kj,K) K 2 (q^ qi ' qj ' K) ) 2 K 2 (q^' K) ) 2 



R s (u) = R s (u), 



(5.39) 



where the definition of R s (u) is taken from (|5.25| ). 
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5.2.1 Summary of k = 3 

The pentagon coefficient is given by 



n+3 

C[Qi, Qj,Qt] = {K 2 f +n Y[ (3 ( f l ' q °' qt ' ,Ps) . 



(5.40) 



s=l 



There are three boxes associated with the triplet of vectors Qi,Qj,Qt- We give the formula for the box 
involving Qi and Qj\ the other two may be obtained by exchanging indices. 



C[Qi, Qj]k= 



(K 



2\2+n 



( 



]ir=i ( P (Qj ,Qi);l I ^ ( u ) I P {Qj ,Q<);2 



n+2 



n+3 



P {Qi,Qi)\l\ K \ P {Qj,Qi)-; 



+ { P (Qi,Qi);i ^ p (0j,Qi);2} > ( 5 - 41 ) 



where we have made the following definitions: 
i? r («) 

Qt{u) 



^ q ^, qj ,K) _ 1){ _ g (c ll , qj ,K) ) + = Q) 



(9& 



qt ■ q^' K) {aM(u) _ + g t ( u = 0), 



and 



(qi,qj,qr,Ps) _ 
s — 

K?e(P s , if, tf t ) + Kh{K h P s , K, K t ) + K 2 e(Ki, Kj, P s , K t ) + ^ 2 e(^, K„ K, P s 



K 2 e(Ki, Kj, K, K t/ 



5.3 The case fc > 4 



(5.42) 
(5.43) 



(5.44) 



If there at least four vectors Ki , then we can use four of them as a basis to expand the momentum vector 
K: 



(5.45) 



i=i 



Using the expression ( 2.1(f ) for Qi(u), we find 



(1,2,3,4) , s 



i^ 2 



K. 



-22- 



Use ( 2.22 ) to define the two null momenta Pi(u) = Q2{u) + XiQ\{u). Then, we have 



(1,2,3,4) 



+ 



41.2,3,4) (Pi(n)|4 1 ' 2 ' 3 ' 4) Q 4 (n) + a { ^ 3A) Q 3 (u)\P 2 (u) 



(Pi(u)|Q 3 («)l*M«)] (Pi(u)\Q 4 (u)\P 2 (u)} (Pi(u)|Q 3 (u)|P 2 (u)] (P 1 (u)|Q 4 (n)|P 2 (n)] 

{Pl(u)\Y.U4^ A) Qi(u)\P2(u)_ 



(Pi («) |Q 3 H I P 2 («)] (Pi («) I Q 4 («) I P 2 («)] 



4 (1,2,3,4) ^2 



Ei=l o. 



(Pl(«)|^|P 2 («)] 



K 2 



Therefore, we can derive the following identity: 

1 

(Pi (u) I Q 3 H I P 2 (u)] (Pi (u) I Q 4 (u) I P 2 («)] 



(Pi(«)|Q 3 («)|P 2 N] (Pi(«)|Q 4 («)|fb(«)] ' 

(5.46) 



K 2 



1 



(1,2,3,4) 



(1,2,3,4) 



+ 



Y:t=i^' 2 ' 3A) Kf-KHPi(u)\K\P 2 (u)] \ v (P 1 ( U )|Q 3 (n)|P 2 (n)] (P 1 (u)\Q i (u)\P 2 (u)} 

Generalizing to our case with k > 4, we have 
1 



UU,i#j( p iiM\QtM\PjiM] 

k 1 

E 



(5.47) 



*=WU (^,i(n)|^|P, li2 (n)] fc - 3 (P^iWIOtWlP^Cti)] s= ^ t Kj^j^a&^K* - K\ 
Thus, we can write 

rfc+ra 



n 



(5.48) 



e n 



a W,t,») K 2 



Urtl (Pr^)\Rr(u)\P 3h2 (u) 



By this formula, we reduce the case with a given n, and > 4, to the case with kgg = 3, and n e g = n+k — 3. 
Using the results of the cases k = 2 and = 3, we see explicitly that indeed the coefficients of boxes are 
polynomials in u. 

The coefficients in ( |5.48| ) can be given more explicitly using ( |2.27| ). We denote the quantity in paren- 
theses by i/sy( K *' K i' K, " Kt ^ anc [ 

U,j,t,s) ^2 is2 



7. 



(Ki,Kj;K B ,K t ) 



Ea=i,i, 



(5.49) 



Kfe(K, Kj,K s ,K t ) + K]t(K u K, K s ,K t ) + K 2 e(K u K,,K, K t ) + K 2 e(K h Kj,K a , K) - K 2 e(K h Kj,K s , K t ) 

K 2 e(Ki, Kj,K, Kt) 



The numerator of ( 5.49| ) is symmetric in Ki, Kj, K s ,Kt; the denominator singles out K s , which is why we 
use the subscript s. 
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5.3.1 The total box coefficient 



We have just shown that when k > 4, we can use ( 5.48| ) to reduce to terms with k = 3, read out box 
coefficients for each of these terms, and eventually add them all up. This approach was useful to prove the 
polynomial property. For computing amplitudes, we would like to carry out the summation once and for 
all. 

We have used ( 5.46D and ( 5.47 ) to derive fl5.48 ). In each term on the right-hand side of (|5.4S| ), there 
is one pentagon coefficient and one box coefficient. The pentagon coefficients are uniquely associated to 
different pentagons, with the various Qts along with Qi and Qj, but the box coefficients all contribute to 
the same box, with only Qi and Qj, so we must add them up. 

Upon inspecting the final expression (5.41) for box coefficients in the case k = 3, we see that our task 
is to check that equation ( |5.46 ) still holds if we replace Qt(u) by Qt(u) and Pi{u) by Pj = Pi(u = 0). Let 
us try to copy the derivation. 



(1,2,3,4) 



+ 



(1,2,3,4) 


(Pi\a^ 4) Q,(u) 


. (1,2,3,4)?; , 
+ «3 V3( u 


)\P2 


(Pi\Qa{u)\P 2 


(p x \Qz{u)\P 2 


(Pl|Q 4 («)|P2 





Pl\Q3(u)\P 2 

Pi\ E-=3«! 1,2 ' 3 * 4) ^(^ = 0) + (a<**)(u) - l)(-go)Et 3 «f ' W) ^|P 2 



Pi\Qz{u)\F2 (Pi\Qa(u)\P 2 



(5.50) 



Then second term in the numerator of ( |5,50| ) is zero because qi ■ go = Qj • Qo = K ■ go = 0. Further, we can 
extend the sum in the numerator of the first term to include i = 1,2, since these terms are individually 
zero when contracted between spinors for P\ and P 2 . 
Because 

4 4 



J2^' 2 ' 3 ' 4) Q i (n = 0) = J2^' 2 ' 3 ' 4) ^- 



we have 



(1,2,3,4) 



(1,2,3,4) 



+ 



D I v^4 (1,2,3,4) rr> D 



Pi\Qz{u)\P 2 (Pi\Qa{u)\P 2 



Pi\Qz{u)\P 2 (Pi\Qa(u)\P 2 



It is now clear that ( 5.46Q and ( 5.47 ) will still hold, if we replace Pji(u) — ► Pji(u = 0) and Qt(u) — > Qt(u). 
Now we perform the summation, using ( PTSP . There are two terms. The first terms are collected to 

give 



(5.51) 





Pji,l(u = 0)\Rr(u)\P ji)2 (u 


= 0) 




(P jitl {u = 0)\K\P ji)2 (u = 


o)] n+2 nti,^(^,i(^ 


= 0)|g t («)|P ii)2 (u = 0) 



The total from the second term from (5.41) cannot be simplified further, since each ptew*J depends 
on both P s and Kf. 
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5.3.2 Results for k > 4 

The box coefficients are given by 



C[Qi, Qj]k>4 



(K 



2\2+n 



11^=1 (• p (Qi,Qi);il- R »'WI^ > (Qj.Qi);2 



^(Qj.QOiil-^I^CQj.OOsa 

-r-rn+A: ,Qf,Ps) 
lls=l P s 



n+2 



Ut=l,t^i,j \ P (Qj,Qi);l\Qt{ u )\ P (Q j ,Qi);2 



P (Qj,Qi);l\ K \ P (Qj,Qi);2 



(Ki,Kj\K w ,Kt) 



P (Qj,Qi);l\Qt( u )\P(Qj,Qi);2 



+{ P (Q j ,Q l );l ^ P (Qj,Qi);2}> 



(5.52) 



(qi,qj,qr,ps) (Ki,Kj;K w ,Ki) 



are denned in ( ^42|) , fljuD, (^Ij ) and ( ^49|) , respectively. 



where _R r (u) , (-u) , f3, 

All u-dependence is inside R(u) and Q(u). This form makes it easier to take the derivative in ( 2.20| ). 
Pentagon coefficients are given by 

T-in+k R (.qi,qj,qr,Ps) 

2\3+n lls=l P> 



C[Qi,Qj,Q t ] = (K 



T-fk (Ki,Kj-,K w ,K t ) ' 



(5.53) 



5.4 The degree of the polynomial 

In the cases k = 2 and k = 3, we have seen explicitly that the maximum degree of the polynomial is 
[(n + 2)/2]. For k > 4, the logic discussed after ( |5.48| ) implies only that the degree is no greater than 
[(n + fc-l)/2]. 

However, we can make a stronger claim by performing a different reduction. In the previous subsection, 
we chose a reduction with a symmetric treatment of the factors (—2£-Pi) in the numerator of ( p.6[ ). 
Alternatively, we could choose not to respect this symmetry in reducing the cases with k > 3. For 
example, we can expand P s in terms of -fQ, Kj,Kt,K and use 

(P Jhl (u)\R s (u)\P Jh2 (u)} = ai q ^' qr ' pi) {P jhl (u)\Q t (u)\P Jh2 {u)} + pj****') (P JlA (u)\K\P Jh2 (u)}{5.54) 

This reduction translates to the following relation: 

Box[m, k, n] — > Box[m — 1, A; — 1, n] + Box[m — 1, k, n — 1]. (5.55) 

Upon iteration, we arrive at /c = 2 for a fixed n. Then, the degree of the polynomial is [(m — k + 2)/2] = 
[(n + 2)/2], which is what we wanted to show. (It is clear that the second term in the above reduction will 
not have a higher degree than the first term.) 



6. Gluon example: A(l+, 2+ 3+ 4+ 5+) 

We now present an application to the five-gluon one-loop amplitude in Yang-Mills theory, which was first 
presented to all orders in e in |32|. By supersymmetry arguments, the computation is equivalent to one 
with a scalar field circulating in the loop. 
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This configuration is totally symmetric, so we need only consider any single, representative cut, say 
C\2- The others can be obtained by permuting labels. The cut integrand within (|2.6| ) is the product of 



two tree amplitudes given in [34], with a factor of 2 for the two internal helicity choices: 



Ji2 = 2A L (-e 2 ,l + ,2 + ,-e 1 )A R (£ 1 ,3 + A + ,5 + j2) 

/» 2 [1 2] V[5|fc 34 li|3] 

(1 2) ((h - k 2 y - n 2 ) (3 4) (4 5) ((4 + k 3 ) 2 - ^)((£ 2 + fc 5 ) 2 - ^ 2 ) 

_ 2 (/. 2 ) 2 [1 2] [5|fc 34 ^3] 

(1 2) (3 4) (4 5) ((£ _ fc 2 )2 _ M 2 )( ^ + fc3 )2 _ M 2 )( ^ + fc34 )2 _ ^2) 

2 (^ 2 ) 2 [1 2] J-2?-i\) 

(1 2) (3 4) (4 5) ((£ _ K x )2 _ M 2 )( ^ _ K2 )2 _ M 2 )( ^ _ K3 )2 _ ^2) • 

where we have defined 

£ = e 1 = k i3 -e 2 , (6.i) 

and also that 

K = k 12 , K 1 = k 2 , K 2 = -k 3 , K 3 = -k u , P 1 = [5 3]k 3 + [5 4]A 4 A 3 . 

According to our definitions in (2.6) and fl2.10| ), we have 

k = 3, m = l, n = -2. (6.2) 

Since triangles appear only when n > — 1 and bubbles appear only when n > 0, we are left with only 
pentagon and box terms in this case. Because m = 1, the degree of the polynomial in u is 0. Thus, we can 
set u = in our formula from the beginning. 

From our formula ( |5.40| ), we find that the pentagon coefficient is 

_ 2(/X 2 ) 2 [l 2]S1 2 a ( qi ,q 2 ,q 3m ) ,„„s 
^pen- (1 2 ) ( 3 4} (4 5p ' ^ 

This coefficient is proportional to (/i 2 ) 2 . The integral I^[(^ 2 ) 2 ] is 0(e), which is what we expect of 
pentagons. 

There are three boxes involved in this cut. Using the formula ( |5.41| ), we find that the box coefficient 
associated with Qi,Q 2 is 

C[Qi Q 2 ] = -(Z^ 2 )^ 1 2 1 / (- P (Q2,Qi);il fl il P (Q 2 ,Qi);2] ^(gimmm) ( P (Q2,Qi)-A K \ P {Q2,Qi);2\ \ 



(1 2) (3 4) (4 5) \ (P( Q2 , Ql y,i\Q3\P(Q 2 , Ql y,2] H1 ( p (QM;i\Qs\ p (QM;2] 
+{ jP (Q2,Qi);l ^ ^(Q 2 ,Qi);2} 

-2(^ 2 ) 2 [1 2] 
(1 2) (3 4) (4 5) 3 



(6.4) 
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in which we have used ( 5.2§| ) to simplify the expression. 



Similar calculations give 



™ r,^ -2(^ 2 ) 2 [12] 

C[QiM = (1 2)(3 4)(4 5) «3 , (6-5) 

r<\n n l — — 2(/i 2 ) 2 [l 2] ( g2 ,g3,gi;pi) / R R \ 

c[g 2 ,Q 3 ] - (1 2)(3 4)(4 5) a 3 • (6-6) 



Now we check these coefficients against the result in the literature [ 34 , |15|] , which is (the factor i / (4ir) 



\2-e 



is omitted, and we also changed the result to our convention of the basis definition (|2.29|)) 



A{1+, 2+ 3+, 4+ 5+) = /1oV /00 w ^ e Lv T^Ts f " 4 *( 4 " 2e ) e (^2, h, k 4 )I° +6 [l] + S2 3S34/f +4,(1) [1] 



e ( 1 ~ e ) / ^-^ 7„_ u. z,^+6 m , „..„. .^+4,(1) 

(12) (23) (34) (45) (51) 

+ S 34S45/f +4 ' (2) [l] + S 45*5l/f +4 ' (3) [l] + S 51 S 12 I? + ^ [1] + -Wf^l] 

- (1 2) (2 3) (3 ' 4) (4 5) (5 1} (^(fei, **, h, h)I^ 6 ] - s 23 s u t (1) lA 

-S 3 4S4 5 /f ' (2) [fl 4 ] - S 45 S 5 ilf ' (3) [//] - S5lSl2^f' (4) [/^ 4 ] ~ S 12 S 23 I®'^ [^} ) . (6.7) 



We have used if [,u 4 ] = -e(l - e)/f +4 [1] and /f [//] = -e(l - e)(2 - e)/f +6 [1]. In Appendix [A], we discuss 
various recursive relations and dimensional shift identities. We now apply the identity (|A.9|) to get 



= (~) + it: (-^) (6.8) 

i=l ' 



Then we see that we should be able to reproduce the following correspondences. For C[Qi, Q 2 ] 

(1 2) (3 4) (4 5) 
for the pentagon, 



(91,92,93^1) = 1 f Sl2S23 + 4ie(kx,k 2 , k 3 , k 4 )-^] ; (6.9) 

3 (1 2) (2 3) (3 4) (4 5) (5 1) V A 5 / ' v ; 



2[1 ^(91,92,93^1) = 1 ( gfc^.fe.As.^J- ) . (G.lOi 



(1 2) (3 4) (4 5) ^ (1 2) (2 3) (3 4) (4 5) (5 1) V As 

We have checked that these equations (as well as the ones derived from the other two boxes) are consistent 
with our definitions. 

7. Discussion 

From the u-dependent formulas for 4-dimensional integral coefficients given in |27|], we have now given 
simpler versions, explicit proofs that they are polynomials, and the formula ( 2.2C| ) needed for the final 



evaluation in D dimensions. In this section, we make some remarks comparing our results with two recent 
papers, pf and [pi. 
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The authors of |30| discussed the calculation of rational terms. The rational contribution may be split 
into two parts (eq. (3) of fjO|l ), namely a term depending on q 2 (which is in our notation) and the 
4-dimensional part. For the former term, the authors of (3^] reduce the calculation into effective Feynman 
diagrams. In our approach, we do not distinguish these two terms; they are treated on the same footing 
by using dimensionally shifted master integrals. To deal with the 4-dimensional terms independent of q 2 , 
proposed the mass shifted method (eq. (16) of p0|| ) and following expansion (eqs. (17), (18), (20) of 
). In our terminology, these are the coefficients of u a , which we have discussed. The proposal of |3(| 
is to choose different values of q 2 , while here we use the derivative. We could also choose to substitute 
numerical values of u, and then find the coefficient from a linear equation, as detailed recently in [46|. In 
another recent paper |5(|, this same numerical approach is implemented in the context of [fjcj ]. 

The authors of [29] treat s 2 (which is u in our notation) as an effective dimension. Thus they are able 
to use a 5-dimensional cut to read off the pentagon coefficient. To get the coefficients, they work in two 
different dimensions, D\ and The paper [^9| has thus given a way to deal with the problem of the 
polarization tensor of a gluon or fermion in arbitrary dimension D. By choosing appropriate loop momenta 
and solving a linear system of equations, they can separate the coefficients into spurious terms and the 
terms with various powers of s 2 . As in our approach, the terms with non-zero powers of s 2 will contribute 
to the coefficients of dimensionally shifted master integrals. 

The methods of [29, [3(1 have been implemented numerically [29, [5(], [5l], |5^[ and have been shown to 



be stable and efficient. 6 We have not yet attempted a numerical implementation of the procedure given 
in this paper, and we leave its assessment to future work. Analytically, our algebraic expressions are the 
most general, since we have not assumed renormalizability, and the power of q in the numerator can be 
arbitrarily high. 
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A. The scalar integrals and dimensional shift identities 

The ©-dimensional scalar integral is defined to be 

dPp 1 



nil ~ [ ] J (2n)Vp 2 (p-K 1 ) 2 (p-K 1 -K 2 ) 2 ...(p + K n ) 2 - l ' j 



6 Note added in revised version: There has recently appeared a numerical implementation of another technique for getting 
rational parts of one- loop amplitudes jiH) , based on the generalized-unitarity formalism of [53| combined with an expansion 



in fi 2 . 
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We will use 55 to denote the dimensionality, so that we can specifically set 

D = 4 - 2e. 



(A.2) 



We also define a very useful symmetric matrix, S, as follows: 

( K 2 {K x + K 2 ) 2 ■■■ (K 1 + K 2 + --- if n _i) 2 \ 



S 



* 

* * 



K 2 
^2 



/T 2 



/0 Si Sl2 • • • S n 

* S 2 • ■ ■ S n i 

* * * s n _i 

^* * * * o y 



(A.3) 



If there are explicit powers of /j, 2 in the numerator, we expand in e as follows, 

W) k ] = 4 D+2fe [i] = -e r(fc) /f +2fc [l] +0(6) 



(A.4) 



and deduce that we only need to calculate the coefficient of the 1/e term of In +2k [1] , in order to get the 
rational term. For bubbles and triangles, we need to consider k > 1; for boxes, we need to consider k > 2; 
and for pentagons, we need to consider k > 3. 

We will use two ways to deal with the higher-dimensional scalar integral, mainly following 0]. 

The first way is by calculating the integral directly, using Feynman parametrization: 



J*[l] = (-irr(n-55/2) 



da-i 



da n 



(A.5) 



where 

n 

a ■ S ■ a = a i a j Sij ■ 

i,i=l 

This integral is easy in the cases of bubbles and one-mass or two-mass triangles. However, for three-mass 
triangles, boxes and pentagons, the integral is complicated. 

The second way is by using a recursive relation, which reduces the higher-dimensional scalar integrals 
to lower-dimensional and lower-point scalar integrals 0: 



l 



(n - 1 - 55) A, 



2tf[l]+5>n,i#? ) [l] 



i=l 



where 



i=i 



(A.6) 



(A.7) 



i=l 



If 55 = 15 + 2k - 2, then we have 

1 



D+2fe 



[1] 



A: 



71 — 3 



i=l 



7n,» \ r D+2(fe-l),(i) 

a7 ' 



[i] 



(A.8) 
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Similarly, we can write 



tD 



2\k 



GO 



k — 1 



n— 3 



A, 



rD 



2\fc-l 



i=l 



7n,i 



T D,(i) 



2\fc-l 



(A.9) 



The recursive relations are very convenient in dealing with three-mass triangle and higher point cases. For 
one-mass and two-mass triangles, the matrix S is singular, so these recursive relations are not well defined. 
However, it is possible to recover the results from massless limits of the three-mass triangle. In practice, 
then, we can always use these recursive relations, taking a massless limit at the end in special cases (also 
boxes and pentagons). 

In the following, we give compact recursive formulas for bubble, triangle, box and pentagon, for 
convenient automated evaluation. 



A.l Bubble 

For the bubble, the matrix S defined in ( |A.3j ) becomes 

s-- l -(° K2 ) 

b ~ 2 [k 2 J ' 

so 

2 



s- 1 



72,1 = 72,2 



K 2 ' 



A, 



\h 



4 



and 



a ■ S ■ a = —a\aiK . 
Using the recursive relation ( |A.8| ), we find 



1 



Alternative, we can use the Feynman parametrization to calculate it directly. 

/f+ 2fc [l] = T{-k + e) f d ai da 2 5(1 - a x - a 2 ) {- ai a 2 K 2 ) k ~ e 
Jo 



T (-k + e) T{k 1 + 1 ~ e) \ (-K 2 ) k -< 
v J T(2k + 2-2e) y ' 

T{k + l)(K 2 ) k 1 Q 
T(2fe + 2) K h 



(A.IO) 



(A.11) 



(A.12) 



(A.13) 



(A.14) 
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A. 2 Triangle 



The matrix S is 





(0 Kl 




1 


(a si s 3 \ 


Kl 




si s 2 






1) 




W S 2 0/ 



(A.15) 



Its inverse is 



5" 



sis 2 s 3 



-«2 S 2«3 «1 
,2 



« 2 \ 



S2S3 -S3 S3S1 

\s 1 s 2 S3S1 -sf/ 



(A.16) 



so 



73,1 
A3 



S2(s2-S 3 -Si) S 3 (S3-S1-S 2 ) 
-, 73,2 = - 



S3S1 



SlS 2 



73,3 



- s 2 - S3) 



S2S3 



s l + s 2 + s l ~ 2 ( s l s 2 + S2S3 + S3S1) 
S1S2S3 



(A.17) 



and 



a ■ S ■ a = — (aio^-fTi + a^a^K^ + 0301-^3). 



(A.18) 



One-mass and two-mass triangles can be evaluated by Feynman parametrization. If K^ = 0, then 

^3?2ro fc [l] = — r(— A; + 1 + e) / daida2daz 5(1 — a\ — a 2 — 03) (—a^Kf — a 2 a3-Kf 

Jo 

T(k) {Kl) k - (Kl) k 1 



-2\k-l-e 



T{2k + 1) Kl - Kl e 



+ O(e ). 



(A.19) 



If, in addition, K^ = 0, then 



T D+2k 
'3:1m 



[1] 



T(k) 



(Kl) k ~ 1 - + 0(e°). 



T(2fc + 1) 

For the three-mass triangle, we use the recursive relation ( |A8| ) repeatedly, and obtain 



(A.20) 



T D+2k 



[1] 



fc-e 



L ^ w) [i]+5E 
-> ' i=i 



A, 



73,i \ jD+2(fe-l),(i) 



[1] 



l^ r(l + l) / 1 rD+«,(i) ri1 , n < 



(A.21) 
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where we make use of our previously derived bubble result. The first several cases, written explicitly, are 
/f+ 2 [l] = -l + 0( e °) 
I° +i ll] = -±- e (s 1 + s 2 + s 3 ) + O(e ) 



/? +6 [l] 



1 



360e 



(si + s\ + s| + sis 2 + s 2 s 3 + s 3 sx) + O(e 



ff +8 [l] = " 6720^ H + s 2 + s i + 4*2 + sis 2 2 + s 2 2 s 3 + s 2 sl + s|si + s 3 s\ + -s x s 2 s 3 ) + 0(e u ) 



Note that A3 in ( A,21| ) has cancelled out of the numerator and denominator. 

We have verified that the results for one-mass and two-mass triangles are consistent with the massless 
limit of the three- mass triangle result. 



A.3 Box 

The matrix S is 



S 



(K 1 + K 2 f K\ 
\ K 2 (K 2 + K 3 f 



K\ (K 1 + K 2 ) 2 K\ \ 
K 2 (K 2 + K 3 ) 2 









/ si S12 s 4 



si s 2 s 23 
S12 s 2 s 3 

\«4 S 23 S 3 J 



(A.22) 



Using the recursive relation (|A.8|) repeatedly, we find 



rD+2fc 



1] 



k-l-e 



A, 



r D+2(fc-l) 



r D+2£,(i) 



7V\ / D+2(fc-l),(i)^ 
A4 



[l]+0(6°). 



Here we use the identities for triangles. The first few cases, listed explicitly, are 



(A.23) 



Jf+ 6 [1] 



1 



rD+8 



[1] 



120e 
1 



2520e 



( Sl + s 2 + s 3 + s 4 + S12 + s 23 ) + 0(e C 



s l + s 2 + s 3 + s 4 + s 12 + s 23 + s l s 2 + S 2 S 3 + S3S4 + S4S1 



+ (S12 + s 23 )(s 1 + s 2 + s 3 + s 4 ) + ^(sis 3 + s 2 s 4 + S12S23) ) + O(e ) 



Notice that the factor A4 in ( A.23j ) cancels out 
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A. 4 Pentagon 

The matrix S for the pentagon is 



S 



Using the recursive relation (I 



1 



k-l-e 

fc-l 



/ Sl Sl 2 S45 S 5 \ 

si s 2 s 2 3 s 5 i 

Sl2 S 2 S3 S34 

•S45 523 S3 S 4 

V s 5 S51 S34 S 4 / 

repeatedly, we find 

5 



(A.24) 



-1 \ T D+2(k-l) 



i=i 



~75/, 



r D+2(fc-l),(i) 
'5-1 



[1] 



^E 



r(/) 



fc-^ 5 



2^r(fe) Va 5 



The first few identities we derive this way are 



D+2t,(i) 



[l] + 0(e°). 



(A.25) 



8=1 



J 5 



[1] 
[1] 



1 



720e 
1 



20160e 



(si + s 2 + s 3 + s 4 + s 5 + S12 + S23 + s 34 + S45 + S51) + £>(e°) 



s l + s 12 + s l s 2 + S12S34 + Sl(si2 + S 2 3 + S45 + S51) + -(S1S3 + S1S34 + S12S23) 



+ four cyclic J + O(e ) 
The factor A5 in ( |A.25| ) cancels out. 

B. Explicit expressions for triangle coefficients 

In this appendix we collect some simplified expressions for triangle coefficients with n < 2. This is sufficient 



in renormalizable theories. For general n, we can always go back to the general formula (3.11). 
For n = — 1: 

rfc-l 



For n = 0: 

a s (K 2 ) 2 



1 / I\j=l (P(q s ,K);l\Pj\P(q s ,K)-fi 



2 \Ut=i^ s (P( qa ,Ky,i\K t \P iqs , K y,2\ 



+ {P(q s ,K);l " P(qs,K);2} ■ 



(B.l) 



C[Q S , K] n= Q 



m=i (p ( 



( ?s ,K);l|- P il- P (g S) if); 2 J 



E 



2pj • <? s 



Afe,*Q \ Uti,t^(P(q s ,Ky,i\Kt\P iqs ,Ky,2] \ U {P^M \ p {<u#)fl] 
2<ft ■ g s 



^);2j 



+ {P(g a ,^);l Pfe,i^);2} 



(B.2) 
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For n = 1, we have linear u-dependence: 

(K*)W S ( X%Z\{P {qt ,K)t\Pj\P( q .,K)V 



C[Q S ,K] 



n=l 



{T\ + ^ 2 ) + {P (gs ^ ); i ~ P (gs ^); 2 } (,B.3) 



where 

^2 



fc+1 

E 



+ E 



(2% • q s 



{ { P {qs,K);l\Pj\P{q s ,K);2\ t=1 ^ s ( P (q s ,K);1 1 K t I P( ?s ,K);2 

2 fc+1 



(B.4) 



fc+i 

^{(P^KyMPjlP^K)?] 



{2pj ■ q s ) 



V fc+1 i / 



\- P (g„«');2|-Pj|-P(g„A-);lJ 



+ E 

t=l,t^s 



<- P (g s ,^);ll^t|^(g s ,^);2] 



t=l,t^8 

Finally, for n = 2, we again have linear n-dependence: 

.2,6 / U-tKP^llPjlP^K)^ 



<x 2 s K* J (P^KyMP.n^] 



C[Q s ,K] n=2 



2 (if 2 

'3A(, s ,K) a Vn t fc =l^ s <^ s ^);ll^l^^);2 



-{Tl + 3^ 2 + ^ 3 ) + {P( q „K);l <- ^.,K);2} , 



where 



fc+2 



(2Pi • q s j 



{2qt ■ q s 



{ {P(q s ,K)-l\Pj\P(q s ,K)-2] t J^^ g {P(q s ,K);l\ K t\P{q s ,K);2\ ) ' 



(B.6) 



fc+2 



-F 2 



E 



(2pj • g a ) 



2 fc+2 



1 \<-P(5.,A-);l|-Pj|-P( fc) A-);2] 

(2g* • ?s) 



gs , A ( P fa^);2|^l^fe,g);l] 



fc 

E 

t=l,t^s 



s P(q„K);l\Kt\P{q„K)fi 



V V" 1 (h ,A ^ I 1^1 ( P (ls,K);2\ K t\P(qs,Ky,l\ ( , 



and 



^3 



fc+2 



E 2 



(2pj • 9 s y 



J 

fc+2 



1 \(-P(g„A-);l|-Pj|-P(g„A-);2] 



E 2 



(2<?t • g s 



t=l,t^s 



(P{q„K);l\Kt\P(q„K);2] 



+ 3 f n _ ,A q * + ^ < P fe^);2l^l p fe^);i] (2Pj-g*) 

+ fen (1 - tt) ^ +1 ; <p (?s ,, );1 |p,ip (?s ,, ); 2] 2 



(B.8) 
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C. Proof of the polynomial property of bubble coefficients 



Here we present a proof that the bubble coefficients are polynomials in u. For this proof, we make use of 



their derivation from spinor integrals [27], along with certain results of Ossola, Papadopoulos and Pittau 
(OPP) [T^] to analyze the integrand. 

Given the cut integral bubble coefficients are given by the sum of residues at the poles of the 

following function (see Appendix B of (27|). 7 

(C.l) 



^ 1 dW n>n „ q (s) 



q\ dsi 



where the residues of B n ^{s) are taken before the derivative in s, and the function B nt t(s) is defined to be 

• w (f|/fif +, (f|.,Aif)"n*,iMep(/f-^)if) 



Bubble contributions in the cut integral ( |2.6| ) appear only if n > 0. In the case that n = 0, it is easy to see 
that we can simply set s = 0. Then, (£\Q p K\£) = -y/T=u(£\q p K\£), and (£\RjK\£) = -y/T^u (£\pjK\£) , 
so the u-dependent factor \f\ — u cancels out of the numerator and denominator. 

Our strategy is to decompose B n: t{s) as a power series in \/l — u. Then we will show that the terms 
with odd powers of \J\ — u correspond to spurious terms discussed by OPP. Therefore, they will vanish 
upon integration, and we will be left with only even powers of \/l — u, i.e. a polynomial in u. 

Note that when we apply the OPP results, we are dealing only with the four-dimensional momentum 
q (or I in our notation), so we do not need any parts of the OPP formulas involving the extra-dimensional 
variable q 2 . Also, in our case we have po = 0, and the mass mf should be shifted to mf + [i 2 . 

We emphasize one point which is crucial for our proof: the one-to-one correspondence between the 
form (2.i ) and the form ( |CJ[ ) in D = 4 dimensions. That is, every factor —21 ■ Pj in ( |2.6| ) corresponds to 
a factor ^|i?j 4Z ^(i^ — sr/)\i^ in (|C.l|) , and vice versa. It is very important that since now we are in pure 



4D, the Rf D) = -Pj, i.e., Rf D) = Rj{u = 0). Similarly for the factor (£\Q p (u = 0)(K - srj)\£) in (|CA|) 

^12 



and factor (p — Kj) in ( p.6[) . In the following proof, we go back and forth freely between these two forms. 
C.l Reducing the number of propagators 

The spurious terms of OPP have at most four propagator factors in the denominator. In order to make 
use of their results, we must begin by reducing our (arbitrary) number of propagators to four or fewer. In 
our formalism, the corresponding condition on ( C.2[ ) is that k < 2, because we have a unitarity cut. 



We perform the reduction in (at most) two steps: first, from k > 4 to k < 3, and then from k = 3 to k < 2. 



7 In this discussion, we drop prefactors independent of loop momentum, as well as the possible prefactor c(/x 2 ). Also, we 
have used (K — srf) instead of (K + srf) in (|4.l|). This change is compensated by dropping the factor (— l) 9 . 
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Reducing from k > 4 to k < 3: 

If k > 4, then there are at least 4 Qj's in the denominator and at least one R in the numerator of 
( |C.2p . Therefore, we can expand the vector R in the basis of the Qi, as follows. 

R = Y j x l Q l (C.3) 

i 

With this expansion, the original term can be expressed as a sum of four others, in each of which there 
is a cancellation between numerator and denominator, reducing k to at most 3. Of course, we must be 
sure that the coefficients x% are independent of u. To see this, expand Qi and R as in ( |2.15| ) and (2.17) by 



writing Qi = — (yl — u)qi + otiK and R = — (yl — u)p + j3K. Then, we find that ( |C.3| ) becomes 

4 

p = '^2x i q i , s ^x i ai = (3 (C.4) 

i=l i 

Here it is clear that the solutions Xi are independent of u. Note that the equation p = Ylt=i x iQi has only 
three independent components, because the vectors qi span the 3-dimensional space orthogonal to K. Thus 
we have four equations giving a unique solution of Xj. 

Reducing from k = 3 to k < 2: 

Now we reduce further, from k = 3 to k < 2. Since k = 3 (and we know n > 0, because we are 
discussing bubbles) there is more than one R in the numerator. Taking any one of the R, we expand 



P = tjrK + y\K\ + y 2 K 2 + IJ3K3. 



Then, 



R = -vT^ (p - ^-K^j + (3K 

= Y^Vi\-VT^i(Ki-^-^K\ + aiA + (/? -Y,Vi^i)K 
i=l ^ ^ ' > i=l 

3 3 

i=l i=l 

Substituting this expansion into the numerator of ( |C.2j ), we obtain 

' ^ 3 s(£\ V K\e) 

Uli(i\Qi(K-sr,)\£) UUU, ¥ tm(K-sr,)\£) lP l) Uli (t\Qi(K - sr,)\£) 

The first three terms fall into the case k = 2. The last term, with the factor s (£\rjK\£) in the numerator, 
still has k = 3, but we see from comparison with (C.2) that we have effectively reduced n by one. 8 
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Repeating the reduction on the last term, n times, we arrive at a term with n = 0. As we discussed in the 



paragraph following ( C.2 ), such a term is independent of u. 

Having accomplished the reduction of our proof to the case k < 2, we proceed to apply the results of 
OPP in a case-by-case analysis for k = 0, 1, 2. 

C.2 Case-by-case analysis 

We now analyze each of the cases k = 0,1,2 in turn, rearranging our integrand (p.2j) so that the terms 



with odd powers of \/l — u take the form of the spurious terms of OPP [17], which were proven there to 
vanish upon integration. 

The case k = 0: 

We apply the OPP result directly and make use of their notation. Recall that we always have po = 0. 
Use pig] ), i.e. Rj(u) = -VT^u (Pj - + fyK, and expand Pj as follows (see (2.23) of g7|): 



Pj = Vl K + y n n + y 7 i 7 + y s l 8 . (C.5) 

More concretely, 

K = £ 5 + 24 (C.6) 
K 2 K 2 

£ 7 = X e5 X £6 , £ 8 = \ £6 \ iB , n = h- 24 (C.8) 

Here r] is the same null vector rj, chosen arbitrarily, that we used inside B n ^{s). We see immediately that 

Rj(u) = -(/3j - y n p n - y £7 /3 £7 - y ts /3 t7 )R K + y n R n + ViRi-, + ysRe 8 (C9) 



where Rg 7 = —y/l — u \ l 7 — ^j^-Kj — ^JfK, and the other three vectors are defined similarly from ( 2.17| ). 
After accounting for orthogonality properties, 

n-K = e 7 -K = £ 8 -K = 

we can see specifically how the w-dependence enters: 



R K = -K, R n = -(VT=u)n, Re 7 = -{VT^)i 7 , Ri e = -(VT=u)£ 8 . (CIO) 



Now we use the one-to-one correspondence between form ( |2.6D and form (|C.2 ), that we emphasized at 



the beginning of this section. In the present analysis of a term of type ( |C.2| ) with k = 0, and a numerator 



8 There is another way to see this point. The presence of a term (£\r)K\£) implies that there is a factor (—2£ ■ K) in the 
form (2.6). By the delta-function condition from the 4-dimensional unitarity cut, this factor is equivalent to K 2 , so we have 



reduced n by one. 
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factor with Rt 7 , the corresponding term in fl2,6|) will have the factor —2£ ■ £j. If we expand every Rj 
according to ( P-9|) , then our general term is of the following form: 

(-21 ■ K) SK (-2£- £ 7 ) S7 (-2£- £ s ) Ss (-2£- n) s ". (C.ll) 

We have introduced integers Si to denote the powers. The ^-dependence of such a term is precisely the 
factor yl ~~ u 

We need to show that if sj + sg + s n is odd, then the term is spurious in the sense of OPP ]l7|| . First 



of all, the unitarity cut condition means we can replace {—21 ■ K) — > K 2 , so the value of n is effectively 
reduced by one, and we can ignore that factor for the rest of the proof. If either sj = or sg = 0, then we 
see immediately from |l7]]-(2.29) that the term is spurious. When both S7,sg are nonzero, we use the ex- 



pression [17j]-(2.33) to reduce to the case of [|17|-(2.29). If sj ^ sg, the conclusion is obvious. But if sj = sg, 



then s n is odd, and so, after applying [17]-(2.33), 2i + s n is still an odd power, and we can again conclude 



with |l7]]-(2.29) that the term is spurious. Finally, we must account for the first two terms in [[T7fl-(2.33). 
The first is £- =o 0(A), which is zero by the unitarity cut condition. The second term is ©(g 2 ), which 
is zero since our present analysis is purely four-dimensional, as we remarked at the beginning of this section. 



The case k = 1: 



We continue using the notation of [17], and also its discussion of 3-point like spurious terms. Now we 



use the following expansion involving the single vector Ki\ 

Pj = y K K + Vi Ki + y 3 £ 3 + y 4 £ 4 . (C.12) 
The vectors ^3,^4 are defined as £3 = A^A^, £4 = A£ 2 A^, where £±,£2 are constructed from K,Ki. Then, 

Rj(v) = -v 7 !^ (Pj - ^^K^j + ftK 

= -(Pj - yioti - y 3 Pi 3 - y^Pe 4 )RK + y-iQi + y?,R£ 3 + y^Ru 

Now we substitute this expansion into (£\Rj(u)(K — srj)\£). The term with Qi cancels the a factor in 
the denominator, returning us to the case of k = 0, which we have already addressed. The term with Rk 
reduces n by one, as discussed above. For the remaining two terms, we use £3 ■ K = £4 • K = to write 
Re 3 = — — u)£ 3 and Ri A = — — u)£ 4 , just as in the case of k = 0. 

However, unlike the k = case, when we put Ri 3 and Rn A back in (C.2), the power of \/l — u is not 



always given by the power of Ri 3 and Ri 4 , so we must be careful. Let us consider the separate cases for 
each term in the expansion. 

• (a) If the term contains neither Ri 3 nor Rg 4 , then either Qi effectively reduces k = 1 to k = 0, or Rk 
reduces n to n = in n steps. Either way, we know from previous analysis that the odd powers of 
\/l — u drop out. 
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• (b) If the term contains Ri 3 or Rg 4 , but not both, i.e., 



2£ ■ h) a (-2£-£ 4 ) b 

or — , a,b^0 



{£ - K t f -tf {i _ K .)2 _ ^2 

by (l7]]-(2.20), the contribution is zero. 

• (c) If the term contains both Ri 3 and it^ 4 ,i.e., 

(C - K,f - v? 

then we need to use the first equation of fT7j-(2.15) to reduce the pair. There are three terms on the 
right hand side of the first equation (remembering that the OiTf') does not exist in our case). The 
first two terms reduce n by two (notice that F depends on /x 2 through the mass), and the third term 
reduces k = 1 to k = 0. By this manipulation, we reduce case (c) to either case (a) or case (b). 

The case k = 2: 

We use the same expansion of Rj as in the k = 1 case, and perform a similar analysis. Factors with 
Rk factor effectively reduce n by one. Factors with Q\ reduce k = 2 to k = 1. For Rp 3 and Re 4 , we need 
to use equation ||l7]]-(2.15) to simplify further. Similar to the case k = 1, we have following three cases: 

• (a) If the term contains neither Ri 3 nor Rp 4 , then either Qi reduces k = 2 to k = 1, or Rk reduces 
n to n = in n steps. Either way, we know from previous analysis that the odd powers of y/1 — u 
drop out. 

• (b) If the term contains R^ 3 or Rg 4 , but not both, i.e., 

(-2i -hr {-2i.h) b 

or , a, 6^0 



((£ _ Ki) 2 - vP)((£ - K 3 f - y?) ((£ - Kif - ^)({£ ~ Ktf - m 2 ) 
then we apply the second equation of |17]]-(2.i5) repeatedly until we reach the form ]17|-(2.18). There 
are three terms on the right-hand side of |l7]]-(2.18). The first term will depend on u polynomially 
through the mass, while the third term will reduce k = 2 to k = 1. The second term is the spurious 
term, which gives zero contribution. 

(c) If the term contains both Rg 3 and Re 4 ,i.e., 

(-2l-£ 3 ) a (-2£-£ 4 ) b 



— , a, 6/0 

{{i _ Ki) 2_ fl 2 ){{i _ Kj) 2_ fl 2 ) 

then we apply the first equation of ||17|]-(2.15) to reduce the pair. Then we use the second equation 
of |]l7|-(2.15) and finally reach the form of (2.18). The discussion of this case is parallel to case (b). 



We conclude that the bubble coefficients, as given in ( |4.1[ ), are indeed polynomials in u. Knowing this 
fact, the degree of the polynomial can then be read off from the formulas of Section 4; it is seen to be [n/2]. 
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D. Pentagon double cut 



The double cut of a pentagon integral, defined according to (|2.6D as 

1 



C[I 5 (K;Ki,K 2 ,K 3 )] 



d^p (pi ) 



: 5(p 2 )5((p - Kf) 



is given by the following expression [21]: 
C[h(K;K 1 ,K 2 ,K 3 )] = - 



du u 



-1-6 



u 



(K 2 f 



+ 



+ 



S[Q3 


Q2, Qi 


K] 


In^ 3 


Q2 - V(Qs 


Q2) 2 


Q3Q2 




Q2) 2 - 


Q3Q2 


Qs 


Q2 + ViQs 


Q2) 2 


— q\q\ 


S[Qs 


Qi,Q2 


K] 


In^ 3 


Qi - ViQz 


Qi) 2 


— Q 2 Q 2 




Q1) 2 - 


Q 3 Qi 


Qs 


Qi + V(Q* 


Qi) 2 


— Q 3 Q 2 


S[Q 2 


Qi,Qs 




m Q2 


■ Qi - ViQz 


■Qi) 2 


-QlQ'i\ 




Q1) 2 - 


q\q\ 


Q2 


■Qi + V(Q2 


■Qi) 2 


- QlQl J 



where S[Qi, Qj, Q k ,K] was defined to be 



S[Qi,Qj,Q k ,K] 



T 2 



with 



-8det 



(K-Q k Qi K Qj K 

Qi ' Qk Qi Qi ' Qj 

\ Qj ' Qk Qi ' Qj Qj 



-4det 



( Q\ Qi' Qk Qj ■ Qk 

Qi ' Qk Qi Qi ' Qj 

\ Qj ' Qk Qi ' Qj Qj 



(D.l) 



(D.2) 
(D.3) 
(D.4) 



(D.5) 



(D.6) 



Here we rewrite (D.5) so that the u-dependence becomes transparent. We need to define a few auxiliary 
quantities. In terms of a particular matrix denoted by S, 



( 

K 2 



K 2 




Kl 



Kl 



\ 



Kl (Ki - Kf (K 2 - Kif (K 3 - K, * 



(Ki — K) 2 (K 2 -K) 2 (K 3 -K) 

3-K1 

K 2 (K 2 -K) 2 (K 2 -Ki) 2 (K 3 -K 2 ) 2 
\ K 2 (K 3 - K) 2 (K 3 - Ki) 2 (K 3 - K 2 ) 2 



(D.7) 
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we define 



/ 



A[Ki;K 2 ,K 3 ,K] = -det 



K 2 K 2 \ 

(K 2 -K 3 ) 2 (K 2 -K) 2 {K 2 -K x f 
K 2 (K 3 - K 2 ) 2 (K 3 - K) 2 (K 3 - K x ) 2 



K 2 




B[K U K 2 ,K 3 ,K] 
C[Ki, K 2 , K 3 , K] 



K 2 (K-K 2 ) 2 (K-K 3 ) 2 
\ 1 1 1 

5 

det(5) ( S ~%> 
2det(5). 




1 



(K - Ki] 
1 



Then, 



S[Q 2 ,Q 3 ,Q U K] 



4K 2 A[K r ,K 2 ,K 3 ,K] 



(D.8) 



(D.9) 
(D.10) 

(D.ll) 



uK 2 B[Ki,K 2 ,K 3 ,K]-C[Ki,K 2 ,K 3 ,Ky 

Now it is evident that the numerator of ( p.ll| ) is independent of u, and the denominator is linear in u. 
Furthermore, B[K\, K 2 , K 3 , K] and C[K\, K 2 , K 3 , K] are totally symmetric in their arguments, indicating 
fundamental pentagon nature. The quantity A[Ki; K 2 , K 3 , K] breaks this symmetry for the first argument, 
Ki, indicating that the corresponding propagator is the one that is eliminated in order to show up as part 
of a box coefficient. 
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